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We consider Lorentzian wormholes with a phantom field and chiral matter fields. The chiral 
' fields are described by the non-linear sigma model wtth or without a Skyrme term. When the 

\ gravitational coupling of the chiral fields is increased, the wormhole geometry changes. The single 

throat is replaced by a double throat with a belly inbetween. For a maximal value of the coupling, 
^ ' the radii of both throats reach zero. Then the interior part pinches oflF, leaving a closed universe and 

' two (asymptotically) flat spaces. A stability analysis shows that all wormholes threaded by chiral 

fields inherit the instability of the Ellis wormhole. 

(N 
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Einstein's theory of general relativity describes gravitational phenomena ranging from planetary motion to cosmol- 
ogy. Among the plethora of solutions of the Einstein equations there are also the wormhole solutions [l| . Wormholes 
are a hypothetical topological feature of space-time that correspond to "tunnels" through space-time. In 1935 Q 
\^ ■ the "Einstein-Rosen bridge" was discovered as a feature of the Schwarzschild geometry. It represents a space-time 
\ manifold consisting of two asymptotically flat universes connected by a throat. Wormholes that connect two distant 
regions of our own universe were discussed by Wheeler in the 1950's 0,11]. However, it was shown that this kind of 
wormholes would collapse if one tried to pass them 0-[^. 
, Wormholes which could actually be crossed in both directions, are known as traversable wormholes. Traversable 
CO ' wormholes were considered in (loj when a new kind of matter was coupled to gravity, whose energy-momentum tensor 
would violate all (null, weak and strong) energy conditions. A suitable candidate would be a phantom field, i.e. a 
scalar field with a reversed sign in front of its kinetic term |ll| - [l5|. W ormholes with such phantom fields have been 
studied in a variety of settings, including stars and neutron stars |l6l - [l8| . However, the presence of a phantom field 
is not necessary to obtain wormholes. For example, they occur in theories when instead of Einstei n g ravity higher 
curvature terms are considered, such as /(i?)-theories or Einstein-Gaufi-Bonnet-dilaton theories |19l - l25| . 

In this paper, we investigate wormhole solutions when a phantom field and chiral fields are minimally coupled to 
Einstein gravity. For the chiral fields we first take the non-linear sigma model, and subsequently allow for the presence 
of a higher order term, i.e. a Skyrme term. Our motive is based on the previous observations that the presence of 
non-Ahelian fields can lead to new interesting gravitational phenomena, since "hairy black holes" were discovered in 
the Einstein- Skyrme model (26l . [27j . 

Here we restrict to static spherically symmetric wormhole solutions connecting two asymptotically flat universes. 
The solutions are characterized by three parameters: the coupling constant to gravity a, the radius of the throat of 
the wormhole, and the value of the chiral fields at the throat. In particular, we study the dependence of the wormhole 
solutions on these parameters. Thereby we find that the geometry of the wormhole changes as the gravitational cou- 
pling is increased. Instead of possessing a single minimal radius, the wormhole develops a maximal radius surrounded 
by two minimal radii, i.e., a double throat forms. At a critical coupling, the two throats reach zero size. Then the 
space-time inbetween the two throats forms a closed universe, which pinches off from the two exterior universes. 
Our main focus is on wormholes that are symmetric under the interchange of the two asymptotically flat universes. 
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We consider these solutions as chiral configurations localized in the vicinity of the throat. In the case of a Skyrmion, 
they might give us an idea of the mechanism taken place when such an extended particle is passing the throat. 
Besides the symmetric wormholes we also consider non-symmetric wormholes, by changing the boundary condition of 
the chiral fields at the throat. 

Stability is an essential question for traversable wormholes. Recently, it was shown that wormholes minimally 
coupled to a phantom field possess an unstable mode [28^,^]. However, self-gravitating Skyrmions possess a conserved 
topological charge and are stable [13]; and likewise, black holes with Skyrmionic hair are linearly stable [3l|. Thus 
one may ask, whether this stability exists also for wormholes permeated by chiral fields. After all they do possess a 
conserved topological charge. However, they may also inherit the unstable mode of the Ellis wormhole jllf . To answer 
this question, we study small spherically symmetric perturbations of the wormholes with chiral fields. We impose a 
harmonic time-dependence on the perturbations and find that all wormholes retain the unstable mode of the Ellis 
wormhole. Thus they are all, linearly unstable. 

The outline of our paper is as follows: In section II, we present the action, the Ansatze and the field equations. We 
discuss the main wormhole features in section HI. The numerical results for the symmetric wormholes are presented 
in section IV. Section V is devoted to the stability analysis of these solutions. We report the main results for the 
non-symmetric wormholes in section VI, and give the conclusions in section VII. 

II. ACTION AND FIELD EQUATIONS 
A. Action 

We consider Einstein gravity coupled to a phantom field and ordinary matter fields. The action 



S = 



r r 1 



^ d'^x (1) 



consists of the Einstein-Hilbert action with curvature scalar TZ, Newton's constant G and determinant of the metric 
g, and of the corresponding matter contributions. These are the Lagrangian of the phantom field 0, 



and the non-linear sigma model Lagrangian 



£ph = ^d^cj^d^-cf, , (2) 



Cch^ —TtIL^L^'} (3) 



where = d^UW and k is a coupling constant. The chiral matrix [/ is a function on the space-time manifold taking 
values in the Lie group SU(2). The last term in the action represents the Skyrme term 

>Csk = ^Tr{F^.F^''} (4) 

with Ffj^y — [L^,L^] and e a coupling constant. 

Variation of the action with respect to the metric leads to the Einstein equations 

G^i, = TZf_,^ - -gfj.iylZ = %ttGT^^ (5) 



with stress-energy tensor 



d£,M 

T^i' = g^i^CM - , (6) 



where Cu = J^ph + J^ch + J^sk is the matter Lagrangian. 

B. Ansatze 

For the study of static spherically symmetric wormhole solutions an appropriate choice for the line element would 

be 

ds^ = -A^dt^ + dri^ + R^dn^ , (7) 
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where dfi^ = d6^ + sin^ Odip^ denotes the metric of the unit sphere, while A and R arc functions of 77. Note that the 
coordinate 77 takes positive and negative values, i.e. —00 < rj < 00. The limits rj — > ±00 correspond to two disjoint 
asymptotically flat regions. 

We parametrize the chiral matrix as 

= cosi^ + isini^ e • r , (8) 

with the unit vector field e 

e = (sin^cos(/3, sin^sin^j, COS0) , (9) 
and the vector of Pauli matrices r. The chiral profile function is a function of rj. 



C. Einstein and Matter Field Equations 



Substitution of the above Ansatze into the Einstein equations GJ'^ = SttGT^ yields 
2RR" + - 1 



i?2 



47rG 



^2(^2^/2 ^ 2 sin^ F) sin^ F{2R^F'^ + sin^ F) 



i?2 



-A + 2RR'A' 



AR^ 
4nG i -d)'^ 



n^{R^F'^ _ 2sin2 F) sin^ F{2R^F'^ - sin^ F) 



AR" + A'R' + A"R 



AR 



AttG 



J,'2 „2 7-./2 



k'F" + 



i?2 



sin-^F 



(10) 



(11) 



(12) 



for the tt, r/r/ and 99 components, respectively. (The LpLp and 99 equations are equivalent.) 

The equations for the chiral profile function and the phantom field are obtained from the variation of the action 
with respect to F and 0, respectively. They read 



2 sin^ F 



AF' 



Asm2F 



[AR^^']' = 



The last equation yields 



i?2 



D 

AR^ ' 



2 ^2 R F + sin^ F 
K R H ^ 



(13) 
(14) 

(15) 



where the constant D is related to the scalar charge. Thus, the phantom field can be eliminated from the Einstein 
equations by the substitution (j>'^ = /A^R'^. 

Next we introduce dimensionless quantities via scaling 



?7=-, , (j)^ K(j) , AttG = — . 

K K 



(16) 



This is equivalent to setting k — 1 and 47rG — a. Finally, we rename fj — rj, R — R, (j) — (j) io simplify the notation. 

We observe that adding Eq. ([TT]) to Eq. (|10p and Eq. eliminates the 0'2 term. The resulting equations can be 
cast in the form 



R" = 
A" = 



A - AR'^ - RR'A' 



AR 



2 sin^ F sin^ F 



2R'A' 



a sin FA 



R 



2^2^/2 ^ gij^2 p 



e2i?4 



(17) 

(18) 
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which form together with Eq. ([T^ a system of second order ODEs to be solved numerically. The scalar charge D for 
the solutions can be obtained from Eq. ([TT|). 
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I?2 = AR^ {AR'^ -A + 2RR'A') + A^ I {R^F'^ - 2 sin^ F) R 



sin^F 



The condition D = const is used to monitor the quality of the numerical solutions. 



'2R^F'^ - sin^ F) 



(19) 



III. WORMHOLE GEOMETRY 



A. Geometry of the Throat 



To obtain wormhole solutions, we assume that the function R does not possess any zero. For asymptotically flat 
solutions, R behaves like \ri\ in the asymptotic regions. Consequently, R possesses (at least) one minimum. Suppose 
R'iVo) — for some ryo and R{r]o) = tq. Then it follows from Eq. (IT7|) that 



R"{m) - - ( l-asin'Fo 
ro 



where Fq = Fijio), and we defined 



sin^ Fo 



(1 - a/aci ) , 



(20) 



sin^ Fo 



2 + 



(21) 



Thus, R possesses a minimum at 770 if a < acr and a maximum if a > acr- The first case a < acr corresponds to the 
typical wormhole scenario: a surface of minimal area separates two asymptotically flat regions. In the second case 
the maximum is a local maximum since in the asymptotic regions R = \r]\. This implies that there are (at least) two 
minima of R, one for 77 < 770 and another one for r] > rjQ. The wormhole then has a double throat. For a sequence 
of alternating minima and maxima the spatial hyper-surfaces would also possess two asymptotically flat regions, but 
they would be separated by a sequence of throats, thus a true multi-throat wormhole would arise. 



B. Boundary Conditions 



For the position of the throat - or beyond acr of the maximal area surface, i.e., the equator - we can choose without 
loss of generality 770 = 0. At this extremal surface - the throat or the equator - we then impose the conditions 

i?(0) = ro , R'{0) = , F(0) = Fo . (22) 

The first condition fixes the areal radius of the throat or the equator, and the second is the extremum condition. The 
third condition fixes the value of the chiral profile function at the throat or the equator, which is a free parameter. 
In this study we will first focus on the special case Fq = 7T.7r/2, when the wormhole solutions are symmetric under the 
interchange of the asymptotic regions. Subsequently, we will consider non-symmetric wormhole solutions. 
In the asymptotic regions we impose the following boundary conditions: 

^(77 ~> 00) ^ 1 , F(77 ^ 00) ^ , F(77 -> -00) -> 77,7r . (23) 

The first condition sets the time scale, whereas the second and third conditions result from requiring finite energy 
and topological charge n. 



C. Topological Charge 



In order to identify the topological charge of the solutions, we consider the chiral matrix J7 as a map of spatial slices 
of the wormhole space-time to the group manifold SU{2) ~ S^. Since U takes constant values in both asymptotic 
regions we can contract each of these regions to a point. The spatial slices then become topologically equivalent to a 
three dimensional sphere, where the north and south pole correspond to the asymptotic regions 77 — > +00 and 77 — )■ —00, 
respectively. Thus the chiral matrix can be regarded as a map between two three-spheres, and the topological charge 
is defined as the degree of the map. With the Ansatz for the chiral matrix Eq. ([SJ and the above asymptotic boundary 
conditions for the profile function F with n ^ 1 the topological charge is equal to one. 
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(a) (b) 

Figure 1: Plots of (a) the metric function R and (b) the chiral function sinF for several values of a for non-linear sigma model 
wormholes. 



IV. SYMMETRIC WORMHOLES 



In this section we focus on wormholes that are symmetric under the interchange of the asymptotic regions. For the 
numerical calculations, we re-parametrize the function i? as i?^ = (ry^ -|- r^)h{rf), since the new function h is bounded 
in the full interval -co < rj < oo, in contrast to R. Then, the boundary conditions i?(0) = tq and -R'(O) = translate 
into h{0) — 1 and h'{0) = 0, respectively. 

To obtain wormholes symmetric with respect to the throat (or beyond Ofcr to the maximum), i.e. 77 = 0, we require 
for the chiral field that cosF is an odd function in 77. This yields F{Q) — Fq = 7r/2 when n = 1. 

The mass M of the solutions is obtained from the asymptotic behaviour of the metric function A, 

A^l-^ , (24) 
V 

where the dimensionless mass parameter fj, is related to M hy fi = aM/Mo, with Mq = Attk. 

The ODEs are solved numerically for the given set of boundary conditions and parameters a and tq. The quality 
of the numerical solutions is good, since the variation of the constant D as computed from Eq. (|19l) is typically less 
than 10"^. 



A. Non-linear Sigma Model (NLS) Wormholes 

We first consider wormholes in the absence of the Skyrme term, i.e. we take the limit I/e^ — 0. In this case the 
set of Einstein and matter equations is invariant under the scalings — >■ A77 and R XR. In order to fix the scale we 
choose vq = I. This leaves only a as a free parameter. 

Observe that Eq. ([TS]) reduces to R^A' = const. However, for symmetric solutions A'{0) = 0. Thus this constant 
has to vanish. Consequently, A — 1 is the only solution which satisfies the boundary condition A{r] — cx3) — >■ 1, which 
implies that the mass vanishes. 

We have solved the coupled set of Einstein- matter equations numerically for < a < 1. As examples we show in 



Fig. [T]the functions i?(??) and sinF(7y) for several values of a. In Fig. 1(a) we see that the areal radius R possesses 
a minimum at ry = when a < aa = 1/2. These solutions possess a single throat. For larger values of a, however, 
R possesses a local maximum at rj — and two minima located symmetrically to each side of the local maximum. 
Consequently, these wormholes possess a double throat. As a is increased the minima become more pronounced. Thus 
the size of the two throats shrinks and tends to zero as a tends to the value one. In this limit the locations of the 



throats are at the values ±7r/2. The behaviour of the chiral profile function is demonstrated in Fig. 1(b) where we plot 
sini^ as a function of rj. Note that the function becomes more and more concentrated on the interval — 7r/2 < < /2 
as a approaches one. 

The limit a — >■ 1 is demonstrated in Fig. [5] We conclude that the limiting solution consists of three parts. On the 
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Figure 2: Same as Fig. 1 but for values of a close to the limit a 



a = 0.99000 
a = 0.99900 
a = 0.99990 
a = 0.99999 



2.5 



N 




3 
P 

(a) 



a = 0.001 00 
a = 0.20000 
a = 0.40000 
a = 0.50000 
a = 0.60000 
a = 0.80000 
a = 0.90000 
a = 0.99999 






(b) (c) (d) 

Figure 3: (a) The isometric embedding of NLS model wormholes is shown for several values of q = 0. (b-d) Three dimensional 
view of the isometric embedding for (b) a = 0.001, (c) a — 0.8 and (d) a = 0.99999. 



interval —tt/2 < rj < tt/2 the areal radius and the chiral profile function are 



In this region the metric reads 



R{r]) = COS77 , F{t]) ^ --r] 



ds^ = -df + dif + cos^ rj dQ^ = -df + dilj , 



(25) 
(26) 
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and describes an Einstein universe, A4e = Rt 'x . The chiral matrix J7 is a one-to-one mapping of the three sphere 
to SU{2) ~ 5'"^. It can be easily checked that the system (l25t together with A = 1 is indeed an exact solution of 
the Einstein- matter equations on the interval —tt/2 < rj < tt/2 [32|. Moreover, the scalar charge D vanishes for this 
solution as can be seen from Eq. (fTQ]). This implies that there is no exotic matter present. 
In the outer regions the solution reads 

Rin) - V-^, F{tj) =0 (?? > |), (27) 

Riv) = -V+l, F{rj)^7r < -|). (28) 

In both regions the metric describes the Minkowski space-time, 

ds^ = -dt^ + dR^ + R^dnl , (29) 

and the chiral matrix U is the trivial map U — ±1. Also in these regions the scalar charge vanishes, as it should for 
a vacuum solution. 

The resulting space-time consists of two distinct Minkowski space-times plus one Einstein universe. The origin of 
one of the Minkowski space-times is identified with the north pole of the spatial of the Einstein universe, and the 
origin of the other Minkowski space-time is identified with the south pole of the S'^. 

Note, that at the "glueing points" (i.e. rj = ±7r/2) the derivative of the function R and the profile function F possess 
jump discontinuities which lead to J-function terms in the second derivatives. However, the metric depends on the 
function R^, which has a discontinuity only in the fourth derivative. Hence the Ricci tensor and the Kretschmann 
scalar possess discontinuities at ?7 = ±7r/2, but no 5-function terms [s^. 

The geometry of a spatial hyper-surface of the wormhole space-times is visualized in Fig. |31 Here we show examples 
of the isometric embedding of the equatorial plane 9 — tt /2 for several values of a. The embedding is given by the 
parametric representation 

p{ri) = R{r,) , z{r]) = Vl - R'^ H . (30) 



We observe that the radius R has a single minimum at z = if a < Ocr, corresponding to the waist in Fig. [3(b)1 The 
minimum becomes degenerate if a = a^x and turns to a (local) maximum for larger values of a, corresponding to the 
Close to the limit a — >■ 1 the outer regions are represented by the two plains whereas the inner 
)here, as shown in Fig. |3(d)[ 



belly in Fig. 3(c) 



region forms a sp 



B. Skyrmionic Wormholes 

We now turn to wormholes in the presence of the Skyrme term. In this case the coupling parameter e can be 
included in the scaled quantities by letting 

?7=^, R= — , (f>^K4) , AttG = , D = De , (31) 
en en k 

or, equivalently, by setting e = 1. 

If the Skyrme term is present, the corresponding Einstein-matter equations are no longer scale-invariant. Then rp 
and a can be considered as free parameters. In Fig. |4] we show examples of solutions for rg — 1 for different values 
of a. We observe that in the limit of vanishing a the metric functions approach the massless Ellis solution, that is, 
A{rj) = 1 and R^ijf) = rj^ -\- r"^. Thus, we obtain a Skyrmion in the background of the Ellis wormhole space-time in 
this limit. As a is increased from zero the metric functions deviate from the Ellis solution. From Fig. |4(b)| we see 
that when a exceeds a critical value, the minimum of the areal radius splits into two minima separated by a local 
maximum. Thus again the single throat splits into a double throat and an equator inbetween. As a becomes large 
the metric function A tends to zero at the equator, while the Kretschmann scalar tends to diverge at the equator. 



In Figs. 5(a) and 5(b) plots of the scaled mass /i/a = M/Mq with Mq = 47rK/e and the scaled scalar charge ^/aD 
are presented versus rg for different values of the parameter a. We observe that for a <1 solutions exist for arbitrarily 
large values of rg. In order to identify the limit rp — > oo we introduce the scaled coordinate fj = rj/ro and the scaled 
function R = R/r^. From the Einstein equations p7)) . p8|) and the chiral equation it follows that in the limit 
ro — >■ c» all terms resulting from the Skyrme term are suppressed as I/^q and will vanish. However, this is exactly 
the limit 1/e^ — > discussed in the previous section. Thus, the Skyrmionic wormholes approach the NLS model 
wormholes (after rescaling) when tq becomes very large, provided that a < 1. 
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Figure 4: Plots of (a) the metric function A, (b) the function R, (c) the chiral profile function F and (d) the Kretschmann 
scalar K are presented for Skyrmionic wormholes when vq = 1, a — 0.01, 0.1 and Qcr ~ 1/3, 1.0, 5.0. 



Let us now consider the case a > 1. In this case, NLS wormholes do not exist any more and we expect a different 
scenario. Indeed, we see from Figs. 5(a) and |5(b)"] that Skyrmionic wormholes exist only up to a maximal value ro^max, 
which depends on a. The mass and the scalar charge tend to zero when ro tends to rg 

In Fig. [S] we demonstrate that in the limit tq 
displays the function R when a 



''o.max the solutions become singular. In particular. Fig. 6(a' 
1.5 for several values of tq approaching the maximal value rg, 



I, max- We see that the 

minimum of R approaches zero at the points ±770. In Fig. |6(b)] it is seen that the Kretschmann scalar diverges at ±770 
when ro tends to the maximal value rg.max- 

The geometry of a spatial hyper-surface of the wormhole space-times is visualized in Fig. [71 Here the isometric 
embedding of the equatorial plane 9 — tt/2 for wormholes is presented when rg = 1 for several values of a (as 
examples). We observe that the radius R has a single minimum at z = if a < acr, corresponding to the waist in 
Fig. [7(b)] The minimum becomes degenera te at a — Ucr and turns into a (local) maximum for larger values of a, 
corresponding to the belly presented in Fig. |7(d)[ 

Incidentally, we have also obtained Skyrmionic wormholes with higher topological charge, but not yet performed a 
systematic study of them. 



V. STABILITY ANALYSIS 



The stability of wormholes is crucial for their physical relevance. Recently, it was shown that the Ellis wormhole 
possesses an unstable mode [H, . The unstable mode was missed before because the gauge condition had been 
taken too stringent. Here we do not follow the approach of [1^ [2^, since this would require the metric in closed 
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(a) (b) 

Figure 5; Plots of (a) the scaled mass /i/a and (b) the scaled scalar charge D^/a are shown for Skyrmionic wormholes as 
functions of ro for different values of ct. 




Tl/rp Ti/ro 

(a) (b) 

Figure 6: Plots of (a) the function R and (b) the Kretschmann scalar K for Skyrmionic wormholes when a — 1.5 for different 
values of ro close to the maximal value ro,max- 



form. Instead, we make a careful choice of the gauge condition, which ensures that we do not miss the unstable mode 
present in the Ellis wormhole. In fact we find that all our chiral wormholes are unstable. 

We restrict to a linear stability analysis in the spherically symmetric sector. Our starting point is the spherically 
symmetric metric in the general form 

ds^ = -e^-H^'n)^Kt:V)de + e-'-^'^'^Uif + e^^^'^^dri^ , (32) 

and the time-dependent Ansatze for the chiral matrix and the phantom field 

U = cosF{t,ri) + ismF{t,ri) e- T , resp. ~ (t){t,ri) . (33) 

We derive the Einstein-matter equations for these Ansatze and consider the perturbed metric, scalar field and chiral 
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Figure 7: (a) The isometric embedding of Skyrmionic wormholes is presented when ro = 1 while a — (Elhs), a — 0.1, a — ctcr 
and a = 1.0. (b-d) Three dimensional view of the isometric embedding of the wormholes when ro = 1 while (b) a = 0.1, (c) 
a = Ocr and (d) a = 1.0. 
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profile functions 



Hv,t) 



Hv)- 



Here F and denote the unperturbed chiral and phantom field functions, respectively, whereas the unperturbed 
metric functions obey 



A. 



1, 



Next we expand the Einstein-matter equations up to first order in the small quantities 5u{'q)^ 5\{'q)^ Sa{ri), S(f>{ri) and 
dF{ri). The resulting set of linear ODEs for the perturbations form an eigenvalue problem with eigenvalue w^. If 
is negative the perturbations increase in time, indicating that the solution is unstable. 

We can use the gauge freedom to reduce the number of linear ODEs. Consider the ODE for the perturbation of 
the phantom field 



R 

— {Scj}' + (j)' {25cT + SX + 6i^)) 



p2 



This form suggests the gauge condition [1| 



6v + 6\ + 2(5(7 = 



(34) 



(35) 



Moreover, we argue that, with (1551) . the function (5( 
is negative. Indeed, taking the gauge condition Eq. 



vanishes identically, if there exists an unstable mode, i.e. if 
5]) into account we multiply Eq. p4p by ^0 and integrate over 
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oOjOo). This yields after integration by parts 



A 



dr] 



Since the left-hand-side of this equation vanishes for normalizable S(j), the integral on the right-hand-side also vanishes. 
However, for negative the integrand is positive. Therefore, the integral can only vanish if Scj) is identically zero. 

With S(t) — and = —6X — 2Sa the resulting equations consist of three second order ODEs for the functions Sa, 
SX and SF in addition to the constraint 



R'SX + 2A 



2a- 



F'SF /2sin2 F 



R 



(36) 



In principle, the constraint for the function 5X can be solved in order to reduce the three second order ODEs to only 
two for the corresponding functions Sa and 5F. However, this would introduce factors 1/R' in the ODEs, which 
diverge when R becomes extremal. Therefore, we choose to solve the system of the three ODEs numerically and just 
verify that the constraint Eq. ([55t is satisfied. 

The boundary conditions follow from the requirement that the perturbations have to vanish in the asymptotic 
regions, that is, 

5a{±oo) = 5X{±oo) = SF{±oo) = . (37) 

In addition, in order to ensure that the perturbations are normalizable we impose the condition (t(0) = 1. The 
eigenvalue ui^ has to be adjusted such that the perturbations satisfy the asymptotic boundary conditions [33 ]. 
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Figure 8: Plots of (a) the scaled eigenvalue for the NLS wormholes as a function of a and (b) the scaled eigenvalue for the 
Skyrmionic wormholes as a function of ro for different values of a. The dot in (a) and (b) indicates the eigenvalue of the Ellis 
wormhole. 

The numerical results are demonstrated in Fig. [8l The case of the NLS wormholes is displayed in Fig. |8(a)[ where 
the eigenvalue is shown as a function of a. Here the dot indicates the eigenvalue of the Ellis wormhole~[28j] . We 
observe that is negative and decreases with increasing a. Thus we conclude that all NLS wormholes are unstable. 
In Fig. |8(b)| the stability analysis of the Skyrmionic wormholes is summarized. In particular, the eigenvalue as a 
function of tq for several values of a is shown. As in the case of NLS wormholes, we conclude that is negative for 
all solutions. Hence the Skyrmionic wormholes are unstable as well. 



VI. NON-SYMMETRIC SKYRMIONIC WORMHOLE SOLUTIONS 



The wormholes we considered so far possess the symmetry A{—ri) = ^(77), i?(— 77) Ri^i) and F(— ?/) = tt — F(r]). 
The last property implies F{0) = Tr/2. However, since F{Q) = Fq is a free parameter we can choose Fq ^ tt/2. Such 
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a choice breaks the symmetry with respect to 77 — >■ —rj and leads to non- symmetric wormholes. As an example we 
consider the non-symmetric Skyrmionic wormholes for fixed parameters Fq — 37r/4 and ro = 1 and varying coupling 
parameter a. The areal radius R and the function sinF are shown for several values of a in Figs. [9(a)] and [9(b)| 



respectively. As can be seen from Fig. |9(a)[ the asymmetry of the areal function R is larger for large a. In contrast, 
Fig. |9(b)| shows that the asymmetry of the chiral field is more pronounced for small values of a. However, one should 
keep in mind that a is the coupling strength in the Einstein equations. Thus, for small a the asymmetry of the chiral 
field does not affect the space-time metric very much, leading to an almost symmetric function R. On the other 
hand, for large a the asymmetry of the chiral field has considerable effect on the metric and induces an asymmetry 
in the metric functions. From Fig. |9(b)| we observe that with increasing a the minimum of the areal function turns 
to a maximum and two local minima next to it. This is similar to the symmetric solutions, however the depths of 
the two minima differ and they are not located symmetrically around the maximum. These non-symmetric wormhole 
solutions thus also possess two non-symmetric throats. 

Fig. 9(c) displays the scaled mass parameter ^/a and the scaled scalar charge Dy/a in terms of a. We observe that 
n/a assumes negative values for small a, possesses a maximum at a ~ 0.26, and tends to some finite value for large 
a. In addition. Fig 9(c) displays the eigenvalue of the unstable mode. Since cj^ is negative we conclude that the 
non-symmetric Skyrmionic wormholes (as considered so far) are, also, unstable. 
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VII. CONCLUSION 



We have considered Morris-Thorne wormholes threaded by chiral fields that carry a conserved topological charge. 
We have focussed on static spherically symmetric solutions, which are symmetric under an interchange of the two 
asymptotically flat iinivcrscs with respect to the throat. 

When the chiral fields are described by the NLS model, the coupling parameter a is the single free parameter for 
symmetric wormholes. For small a the chiral fields have hardly any influence upon the wormhole geometry, and the 
solutions possess a single throat. As the coupling parameter increases, the presence of the chiral fields starts to be felt 
by the geometry. At a critical coupling ctcr, the throat becomes degenerate. Beyond acr the wormhole then exhibits 
a set of three extrema, a (local) maximal surface and two minimal surfaces located symmetrically, one on each side. 
Thus the wormhole possesses a double throat with a belly in the interior. 

As a increases further the chiral fields become more and more localized in the inner region of the belly. At the 
maximal coupling amax finally, the chiral fields arc fully localized inside this inner region, while the areal radius of 
the two throats decreases to zero. The space-time then splits into three parts. The outer parts correspond to empty 
Minkowski spaces, while the inner part corresponds to an Einstein universe with a chiral field carrying topological 
charge one. 

In the presence of a Skyrme term, we have one more free parameter. Thus we can vary the gravitational coupling 
and the throat size independently. With increasing a we again observe an increasing infiuence of the matter on the 
geometry of the throat, and the formation of two throats with an inner belly beyond a critical value acr- The splitting 
of the space-time into three parts is, however, less smooth than in the NLS case. 

By choosing the boundary conditions asymmetrically with respect to the throat we can distribute the chiral fields 
asymmetrically with respect to the two asymptotically flat universes. Then the deformation of the wormhole by the 
matter is asymmetric as well. Consequently, when the single throat splits into a double throat at a critical coupling, 
the two throats develop asymmetrically with respect to their location and size. This may possibly lead to only a 
single throat reaching zero, resulting in a breakup of the space-time into only two parts at a maximal coupling. 

Let us now think of the Skyrmionic wormholes as Skyrmions passing a wormhole. Clearly, one would have to make 
time-dependent calculations to model such a transit. However, the static calculations may already give us some hints 
of what to expect. In particular, we conclude that the Skyrmions will deform the wormhole geometry and possibly 
even let the space-time pinch, resulting in two or three disconnected parts. 

Let us conclude by stating that the chiral flelds cannot stabilize the wormhole space-time. The instability of the 
Ellis solution is inherited by the wormholes threaded by chiral fields, although they do carry a topological charge. 
According to our observations this is not surprising though, since the topological charge may finally simply reside in 
a single of several disconnected space-time parts. 
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